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Abstract 



We study a dynamic model of relaxor ferroelectrics based on the spheri- 
cal random-bond — random-field model and the Langevin equations of motion 
written in the representation of eigenstates of the random interaction matrix. 
The solution to these equations is obtained in the long-time limit where the 
system reaches an equilibrium state in the presence of random local electric 
fields. The complex dynamic linear and third-order nonlinear susceptibili- 
ties an d X3( w )) respectively, are calculated as functions of frequency 
and temperature. In analogy with the static case, the dynamic model pre- 
dicts a narrow frequency dependent peak in xs(T,w), which mimics a tran- 
sition into a glass-like state, but a real transition never occurs in the case 
of non-zero random fields. A freezing transition can be described by intro- 
ducing the empirical Vogel-Fulcher (VF) behavior of the relaxation time r 
in the equations of motion, with the VF temperature To playing the role of 
the freezing temperature TV. The scaled third-order nonlinear susceptibility 
a' 3 (T,u>) = X3(cj)/x / 1 (3a;)x / 1 (a;) 3 , where the bar denotes a statistical average 
over To, shows a crossover from paraelectric-like to glass- like behavior in the 
quasistatic regime above TV. The shape of and x^iuj) — and thus of 
a' 3 (T,u>) — depends crucially on the probability distribution of r. It is shown 
that for a linear distribution of VF temperatures To, a' 3 (T, uj) has a peak near 
TV and shows a strong frequency dispersion in the low temperature region. 

PACS numbers:77.84.Dy, 64.70.Pf, 77.22.-d 
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I. INTRODUCTION 



Relaxor ferroelectrics (or relaxors) represent a new low-temperature state of polar di- 
electrics, which can be regarded as an intermediate state between dipolar glasses and normal 
ferroelectrics |L|,||. Some of the concepts developed for dipolar glasses, such as the Edwards- 
Anderson (EA) order parameter, are applicable to relaxors as well, as recently shown for 
PbMg 1/3 Nb 2/4 3 (PMN) [1, PbSc 1/2 Ta 1/2 3 (PST) i, and Pb^La^Zr.Ti^Oa (PLZT) 
0. In contrast to dipolar glasses, where elementary dipole moments exist on the atomic 
scale, the relaxor state is characterized by the presence of nanosized polar clusters of vari- 
able sizes. This picture constitutes the basis of the superparaelectric model H and of the 
more recent reorientable polar cluster model of relaxors |2|||. By including explicitly the 
long-range frustrated intercluster interactions of a spin glass type into this picture, one 
arrives at the so-called spherical random-bond — random-field (SRBRF) model of relaxor 
ferroelectrics ]7[], which is capable of describing the static behavior of relaxors, such as the 
line shape of quadrupole perturbed NMR in PMN || and PST |4[], and the sharp increase 
of the quasistatic third-order nonlinear dielectric constant [0J|. 

The unusually large value of the static linear dielectric permittivity can also be explained 
within the framework of the SRBRF model if one assumes that the mean value of the random 
coupling J is very close to — but slightly smaller than — its r.m.s. variance J, whereas 
in dipolar glasses the latter is usually dominant. By including an ad hoc electric field 
dependence of Jo into the model, one can furthermore describe the transition from the 
relaxor to an inhomogeneous ferroelectric state for fields E exceeding a critical value E c 
||. It should be noted, however, that the random electric fields, which exist both in dipolar 
glasses [10] and in relaxors @J^,§], seem to be much weaker in the latter case. It is interesting 
to note that Korner et. al [|TTJ have reported a relaxor phase in the dipolar glass DRADP 
in a narrow range of concentration just above the dipolar glass phase, where the system is 
expected to behave as an inhomogeneous antiferroelectric. 

While the static SRBRF model describes a relaxor system in thermodynamic equilibrium, 
there are a number of phenomena suggesting that relaxors, in particular their low temper- 
ature state, are dominated by nonequilibrium effects. Typical examples are the difference 
between the field-cooled and zero-field cooled static dielectric constant, and the occurrence 
of strong frequency dispersion in both the linear and nonlinear dielectric permittivity at 
low temperatures. Is is clear that these properties can only be discussed within a dynamic 
model. In the present paper, we introduce a dynamic model, which an extension of the 
SRBRF model to dynamic problems. Following Vugmeister and Rabitz we assume that 
polar clusters can reorient with a characteristic relaxation time r and write down the corre- 
sponding Langevin equations of motion, which are based on the static SRBRF Hamiltonian. 
These equations explicitly contain the frustrated interactions between the polar clusters 
and thus allow us to study the effects of these interactions on both the equilibrium and 
nonequilibrium properties. In particular, we will discuss here the anomalous temperature 
dependence of the nonlinear dielectric susceptibility and the crossover from the paraelectric- 
like to inhomogeneous ferroelectric-like behavior observed in PMN and PLZT |5|]. As in spin 
glasses, the Langevin equations based on the spherical model can be solved exactly [JT2 



some additional features appear in view of the presence of random fields JT3|. Here we will 
focus on the asymptotic solutions corresponding to equilibirum dynamics, such as observed 
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in a typical dielectric relaxation experiment. 

In Sec. II we introduce the uniaxial SRBRF model Hamiltonian in the representation 
of eigenstates of the random interaction and write down the Langevin equations of motion. 
The asymptotic solution is studied in Sec. Ill, where the static linear and nonlinear suscep- 
tibilities are derived. In Sec. IV the dynamic linear response is given, and in Sec. V the 
corresponding results for the third-order nonlinear response are derived. Finally, in Sec. VI 
we present our conclusions. 



II. DYNAMIC SRBRF MODEL 



In general, the polarization of i-th polar cluster, i = 1,2, ...,N, is a three component 
(n = 3) vector Si = (S ix , S iy , S iz ), its length being restricted solely by the spherical condition 
Yli(Si) 2 = 3iV. In the present work we will discuss the simpler uniaxial (n = 1) case 
— vN < Si < +VN, where Si is subject to the spherical condition 



N 



"£S? = N. (1) 



1=1 

The SRBRF model Hamiltonian is thus 

n = JijSiSj - h i& -gEj2 s ^ ( 2 ) 

ij i i 

where are the randomly frustrated intercluster interactions, h{ local random electric 
fields, E an applied uniform electric field, and g the appropriate dipole moment 0. As 
usual, Jij is assumed to be infinitely ranged and distributed according to Gaussian statistics 
with average value Jo/N and cumulant variance J 2 /N. The Gaussian random fields hi are 
characterized by the random average 

{hihj}^ = A5ij . (3) 

The uniaxial SRBRF model (Q) has potential applicability to uniaxial relaxors such as 
Sr 1 _ x Ba x Nb20 6 (SBN). The present results can be, however, generalized to the isotropic 
n = 3 case as long as there is no mixing of the x, y, z components 0. 

The Langevin equations of motion for the variables Si(t) are written as 

T m = -«m- Miim+ M. (4) 

Here r is the characteristic relaxation time for the reorientation of polar clusters. Eq. (|J) 
implies that r is site independent, however, some variation of r across the system should in 
principle not be excluded, resulting in a distribution of relaxation times [0,0 . The function 
z(t) plays the role of a Langrange multiplier enforcing the spherical condition ([!]) at all times 



The stochastic Langevin forces ensure the proper equilibrium distribution and are 
determined by their ensemble averages 
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m)Ut'))av = 2T5 lJ 5(t-t') . (5) 

Following the theory of spherical spin glasses we now transform to the representation of 
eigenstates ip\(k) and eigenvalues J x of the random matrix J^- [p!4[|12| , ^3| . This is done in 
two steps fll5| : First, one introduces "spin wave" states Sk = iV^ 1 / 2 ^ exp(ikRi)Si; next, 
these are expanded in normal modes 

s x = J2Mk)s k . (6) 

k 

The transformed equation of motion ([|) becomes explicitly 
dS x 



dt 



f3[J x - 2z(t)]S x + (3h x S x + gV x (0)E(t)S x + ■ (7) 



Here \I/a(0) = Ni[)\(Q) and we have rescaled the time to a new dimensionless variable t — > t/r. 
Assuming a field E(t) applied at t = and introducing the integrating factor 



(j) X (t) = exp 

we obtain the solution 



(3J x t-2 [ dt'z(t') (8) 
Jo 



Sx(t) = Mt)Sx(0) + / dhpj^r \j3h x + PgVMEfa) + £ A (*i)] . (9) 

Jo <P\{ti) 

The correlation function 

A 

must satisfy the equal time relation C(t, t) — 1 at all times in view of the spherical condition 
(|). From Eqs. (@)-(0) with the aid of Eq. (|) we thus find 



This is an implicit equation for the Lagrange multiplier z(t). The two types of averages are 
defined as 

«A»o = jfY,f* = I dJ *P*WfW ; (12) 



= a? E *A(0) a /A = / P( Ja)/( Ja) , (13) 

A 

where Po{J\) and p( J\) are the densities of eigenvalues in the k ^ and = sector of the 
spectrum, respectively The eigenvalues J x have a continuous spectrum —2J<J X < 2 J. 
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Jo + J 2 / Jo l|14| . Here we will only 



If | Jo | > J, there is also a discrete eigenvalue at J m 
discuss the case | Jq\ < J. The k ^ density of states is given by the Wigner semicircle law 

MM 



Po(J\) 



1 



2n,P 



(4J 2 - Jl) 



2\l/2 



The k = sector, on the other hand, has the density (15 

1 



2tt(J 2 + Jl - J J)A) 



(4J 2 - Jl) 



2\V2 



(14) 



(15) 



This density of states has a statistical weight 0(1/ N) and is thus relevant only in averages 
containing factors of the type \I/a(0) 2 oc O(N). 

The dielectric polarization of the system can be expressed in terms of the solution @ as 

P W = A? E ^A(O)0 A (t)S A (O) + (3g 2 jf dfi (^j) • (16) 



As shown by Cugliandolo and Dean ||13|| , for times larger than a limiting time t c the 
system in which A^O will always reach an equilibrium state and will thus be characterized 
by equilibrium dynamics. All information about the initial state S'a(O) is lost for t ^> t c , 
i.e., the first term in Eq. ( P!B[ ) becomes irrelevant. In the present case, t c is estimated as 
t c pa 2tJT/A. Typically, the asymptotic regime t 3> t c is explored in a dielectric relaxation 
experiment. In the following, we will limit ourselves to this regime. Also, for simplicity we 
will henceforth set g — 1. 



III. STATIC DIELECTRIC RESPONSE 

We first consider the case of a static electric field E(t) = E applied at t — 0. At 
asymptotic times t/r ^> 1 the system reaches equilibrium and the Lagrange multiplier z(t) 
tends to a constant value z. Thus the function (^) becomes 

Mt)~e- gxt l (g x = 2z-(3J x ) , (17) 



and we can evaluate the integrals in Eqs. (|TT|) and (|T6|). Assuming that 2z > f3J\ for all A 
(to be justified later) we derive the equation for z: 

1 ^ «^)) + ^ ( w^l + ^ (w^) • <"> 

The static linear susceptibility xi — (dP/dE)E=o is derived from Eqs. (|TB|) and (|T7|): 

(19) 

The averages in Eqs. (0) and (^) can be expressed in terms of the generalized averages 
obtained by adding an imaginary generating field iy to the variable g\, namely, 
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(g x - iyf n + x ) 



These averages can be evaluated with the aid of Eqs. (|T2"|)-(|H)|) for n ■ 
n times with respect to iy, and setting y = 0. For example, from Eqs. 
n = and y = we find: 



xi = xf(o) 



z — r — pj 



(3(J 2 + J 2 ) - W-hz 



(20) 

0, differentiating 
[SI) and fl20|) with 



(21) 



where r = \/ z 2 — j3 2 J 2 . 

The n = 1 average is given by 



{1] = z(y)-r(y) \ 1 _ 2^k 
Xl [V) 2 D(y) [r(y) D{y) 



+ 



(22) 



where z(y) = z - iy/2, r(y) = y/z(y)* - (3 2 J 2 , and D(y) = fi\.l 2 + J 2 ) - 2(3J z(y). 
The above equation for z, Eq. fll!^) , becomes in this notation: 



xr J (0)o + /3A x S 1] (0) + /3EM iJ (0) 



2..W 



(23) 



where the averages xi°'(0)o an d X^i\fyo are obtained by setting J = in Eqs. (|2l|) and (p2] 
respectively, and Xi (0) is given by Eq. ( |22"D with ?/ = 0. 
We will also need the n = 2 average 



xf ] (0)c 



8/3 J 2 



(24) 



A numerical solution z(T) of Eq. (^) in zero field (i? = 0) can be found at all temper- 
atures and is independent of Jo as long as |Jo| < J. An example is shown on Fig. 1 for 
A/ J 2 = 0.001. The inset shows that z — (3 J is always positive, and since 2 J is the largest 
eigenvalue of J\ one can see that indeed 2z > 2j3J\ for all A as assumed earlier. When both 
E 7^ and Jo 7^ 0, there are in general two complex solutions for z(E,T) and the present 
theory is not applicable. In the following we will only consider the cases in which a real 
solution z(T) exists and has a real second derivative z" = d 2 z/dE 2 at E = 0. 

For A = and E — 0, Eq. (|23|) reduces to the equation derived in Ref. which has 
the solution zq — (1 + (3 2 J 2 )/2 for T > J. For T < J, however, the solution does not exist 
and z must be obtained from the saddle-point condition [TJJ], yielding z = [3 J '. 



A numerical evaluation shows that the expression (21) for the static linear susceptibility 
fully agrees with the result obtained by means of replica theory in Ref. [0. 

One can also calculate the static third-order nonlinear susceptibility X3, which is defined 
in terms of the expansion 



P = XiE- X3E 3 + 
Obviously, %3 = — {^-/6){d 3 P/dE s ) E =o- Using Eqs. 



(25) 



FT)) and © we find 
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X3 = xi 1] (0K'. (26) 
Evaluating z' ' = (d 2 z / dE 2 ) e=q from Eq. ( |23| ) we obtain the result 

v _ #2 ^1 ] (°) 2 /o 7 \ 

X3 xi 11 (0)o + 2^A x ri(0)o' ^ J 

where Xi (0)q is given by Eq. ([24]) above. The last expression has been evaluated numeri- 
cally and found to be precisely equivalent to the result derived in Ref. using the replica 
formalism. 



IV. DYNAMIC RESPONSE 

We now consider the case of an oscillating electric field E(t) = Encos(ut). This is 
inserted into Eq. (|i~6|) . At asymptotic times t^> t c the response can be written by analogy 
with Eq. ([§) as [I§ 

P(t) ~ [P u e~^ + P-^e-^ 1 + •••]+ c.c. , (28) 

where P^ and P^ are the amplitudes of the first and third harmonic response, respectively, 
which are given by 

P» = xi,oM (f ) + xi» (f ) 3 + • • ■ ; (29) 



En 



Pzu = X3,o( u ) ( -y J +X3,i( w ) 




(30) 



Here we have introduced the linear dynamic response Xi,o(^)> the third-order nonlinear 
responses xi^ii 00 ) an d X3,o(^)> etc. We will focus on the first harmonic linear response 
Xi,o(w), which is equivalent to the dynamic linear susceptibility Xi( u ) — Xi,o(u), and on 
the third-order nonlinear response X3,o(w). The latter is typically measured by monitoring 
the third harmonic component of P(t) at small amplitudes of the field Eq In order 

to ensure the proper static limit u — > we will define the third-order nonlinear dynamic 
susceptibility as Xs{ u ) = ~X3,o( UJ )- From Eqs. (|28| ) - (|30|) we thus find 



OP, 



d(E /2) 



(31) 



E =0 



1 <9 3 P, 



3u> 



6<9(£ /2) 3 



(32) 



£ =0 



In the asymptotic regime, the function 4>\(t) in Eq. (|IE| ) behaves as 



<M*) 



~ e 



(33) 
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where we have defined ip(t) = J" * dt'[z(t') — z], with z representing the solution of Eq. (fL8|) . 
The first part of the response (TTOT) , which will be proportional to ~ E exp(— iut), is now 
given by 

P(t) = (3 f^j ^^(e-^^Ole^bW-^tOle-^ + cx . (34) 



A. Linear dynamic susceptibility 

The part of -Pi(cu), which is linear in Eq, is trivially obtained from Eq. fl34|) by noting 
that <f(t) = for £"0 = 0. We can thus evaluate the integral and using Eq. (|3l|) we find 



^'W^W- (35) 

Comparing with Eqs. (|T9|) and ( p0|) we realize that the averages of the above type can be 
evaluated with the aid of Eq. Q2DD, in which we set y = uj and n = 0, yielding (with r 
restored) 

z - - y/(z - IUT/2Y - (3 2 J 2 - (3J 

Xl[Uj) ~ P(J* + J*)-2J (z-iu;T/2) ' (36) 

For uj — > this obviously reduces to the static susceptibility fl2ip . 

The temperature behavior of wm crucially depend on the temperature variation of 

the relaxation time r(T). The SRBRF model (0) and the equations of motion (|j) contain no 
information about r(T). It has been found empirically |l7|p|,|6| that some of the properties 
of relaxors can be described by assuming a Vogel-Fulcher (VF) relationship for r, namely, 

r = t exp ( jT^jU ' ( 37 ) 

where To is the VF temperature. This expression is valid for T > Tq and would lead to a 
divergence of r for T — > T . There is no a priori relation between T and the parameters of 
the SRBRF model. A similar situation occurs in Ising dipolar glasses, where a probability 
distribution of relaxation times g(lnr) has been used in combination with an empirical 



Debye-type response |18[]. With r lying in the range r min < r < r max , the VF temperature 
T has been identified with the freezing temperature Tf. On the other hand, r m j n has been 
fitted to an Arrhenius-type expression r min oc exp(E/T). The same approach was found to 



be applicable to relaxors as well |L9 



An alternative approach is based on the master equation for the reorientation of cluster 
polarization assuming a VF relaxation time of the type (j37f), where the barrier heights U 
are distributed according to a Gaussian probability function @,^]. Such an approach was 
found to be applicable to PMN and PST in the region T > T . 

In general, we can thus introduce the average dynamic susceptibility 



f Tmax dr , . . , , 

/ — (/(iBr)xiH, (38) 

" Train 
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where the probability distribution of relaxation times g(ln r) is physically justified by the fact 
that relaxors are inherently inhomogeneous systems due to compositional disorder. Thus one 
may imagine, for example, that the relaxor system consists of a set of macroscopic regions, 
which are formally characterized by the same microscopic equation of motion, but differ in 
the value of the parameter r. 

One encounters serious difficulties in attempting to describe the dynamic response at T < 
To. Formally one could assume that r — > oo for T < Tq, but this will lead to a zero value of 
Xi(co>) = at all temperatures T < T . We can single out the following representative cases: 
(i) a single VF relaxation time ([IT]); (ii) a nonsingular distribution of barrier heights g(U); 
(iii) a distribution of relaxation times g(lnr) such that its normalization j Tmax dr g (In r)/r 
diverges as r max — > oo. The first case is illustrated in Fig. 2, where we show the calculated 
real and imaginary parts of xi(T,u) = for several values of frequency assuming a single 
VF relaxation time fl3~7l). As in Fig. 1 we assume Jo = 0.9 J and A/ J 2 = 0.001, as well as 
T = J. Such behavior is incompatible with experiments, which generally show a smooth 
decrease of x'i(T,uj) and Xi(T,uj) across the region where T is expected to be located. 

A more realistic description can be obtained, for example, by assuming a distribution 
w(T ) of VF temperatures T , where T is allowed to vary in the range < T < T™ ax . 
Using the relation d(\nr)g(lnr) = dT w(T ) in Eq. fl3"5| ) and choosing a linear distribution 
w(T ) = 2(1 — T / J)/ J with T™ ax = J, we obtain the temperature dependence of the linear 
susceptibility shown in Fig. 3. Here we used the same set of parameters as in Fig. 2. In 
contrast to the single VF temperature case, the above distribution leads to nonzero values of 
Xx{ u ) a t ai l temperatures. The shape of the real and imaginary part of Xii 00 ) is m qualitative 



agreement with the observed relaxation spectra in PMN |19| and PLZT |[20|| . 

It should be noted, however, that the above result for the linear susceptibility contains 
only the contribution of polar clusters. Other contributions may exist — for example, that 
of optic phonons — which are not expected to show any anomalies near Tf. In general, such 
contributions can be written as a sum of Debye-like terms, with the possibility of an average 
over the corresponding relaxation times. At present, the problem of a realistic relaxation 
time distribution in relaxors, which would be appropriate at all temperatures, has not yet 
been resolved. 



B. Zero-field-cooled susceptibility 



In analogy with Ising spin glasses |2T| and dipolar glasses [22| one can introduce an 
effective relaxation time for the low-frequency response 

-i 



-id lim 



dxiiu) 
duj 



(39) 



Returning to Eq. (|3l"l) and using the above definition of Xi{v) from Eq. (|22"D in which we 
set y = ojt, we obtain the result 



'eff 



—3 t lim 



(40) 



It will be shown later that the real part of the nonlinear susceptibility X3{ u ) is a l so pro- 
portional to x\ (^0 an d in the static limit shows a sharp peak at the "freezing" temperature 
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Tf ~ J for A <C J 2 , while it actually diverges if A = 0. Since Xi{ u ) * s a wen behaved 
function of temperature, it follows that r e ff increases as ~ (T — J) -1 on approaching Tf, 
however, it remains finite at Tf. Thus the behavior of T e ff mimics the freezing transition 
in this case. A true freezing transition at Tf = T can be described by assuming a VF 
temperature dependence of r, which is then transferred to r e // via Eq. (|40|) . 

When r e ff is large, the system will only reach equilibrium at times t ^> r e ff, which may 
become very long. In measuring the static susceptibility xi one should distinguish between 
the field-cooled and zero-field-cooled susceptibility, Xi an d Xi FC i respectively. Here X\° 
is given by Eq. (|TUp and corresponds to an experiment carried out on a time scale t ^> T e ff- 
On the other hand, xf FC can be obtained by turning on a field E(t{) = E at time ti = 
and measuring the response P(t) at time t > t\ as described by Eq. (|l~6l). Thus we can write 



x ^ (t) ^^ Jodtl{m} , (41) 

and using Eq. (|TTf) we find 

// 1 _ p-{ 2 *-0J\)t/T\ 

*r°®='{ 1 L-m ))■ < 42 » 

For t — > oo this reduces to the previous result (|19|) , which corresponds to xf C - The difference 
between the two susceptibilities Sxxif) = Xx ~ X\ FC (t) is given by 

// -(2z-/3J x )t/T \\ roo 

^»=K( W)='l • (43) 

The value of the integral can be estimated by first using Eqs. (|13|) and (]15[) and expanding the 



integrand in powers of Jo/ J To leading order the result is independent of the parameter 
Jo and shows that for A C J 2 and T < J one has a power law behavior Sxi(t) ~ (V r ) -1 ^ 2 ) 
implying a slow decay and a large difference between the two susceptibilities, which has 
been observed experimentally |19[. On the other hand, for T > J the asymptotic behavior 



is a combination of power law and exponential, i.e., 5xi(t) ~ (t/r) -3 / 2 exp[— 2(z — (3J)t/r\. 
Thus in this regime the difference decays much faster and the two susceptibilities become 
indistinguishable on a typical experimental time scale. 



V. NONLINEAR DIELECTRIC SUSCEPTIBILITY 



To calculate the third order partial derivative in Eq. ( p2|) we return to Eq. (g4j), in which 
<p(t) is now a function of Eq. In general, (p(t) will be a sum of terms, which are even powers 
of Eoe ±tujt . We will focus on the second-order term ~ EQe~ 2tujt . Introducing the function 



X(t) 



dMt) 



d(E /2) 2 

we can express the third partial derivative in Eq. 

E o =0 



E =0 

2D as 



d{E /2f 



-Qf3e 



2iujt 



dt ie *w(*-*0^ e -ffA(t-ti)^[x(t) -X(h)] . 



(44) 



(45) 
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The function X(t) will be calculated from Eq. (|TT|) in the asymptotic limit. Considering 
only the terms, which asymptotically behave as ~ e~ 2lu3t and taking the second derivative 
with respect to E /2 leads to 



= -4 

+(3 2 



dt 1 {e-**W\-2(PA / dh / dt 2 ( 



dti 



-2g x (2t-t!)-t2 



[X(t)-X(U 



dtAe-^-^-^e-^^ 



The last double integral becomes for asymptotic values of t 

1 



-2iuit 



(g x - iuj) 2 

We now apply the Laplace transform to Eq. (HBf) using the definition 



X(p) 



dte~ pt X(t) 



(46) 



(47) 



(48) 



and obtain the result: 

P 2 



X(p) 



i 



(g\-iu>) 2 



2(p + 2iuj) 



9\ 



+W/o 



+ {3 2 A 



9l //(, 



a\(a\+p) //o 



(49) 



The averages can be expressed in terms of the generalized responses fl22|). 

The function X{t) can be obtained by the inverse Laplace transform. Its behavior is 
determined by the poles p k = p' k + ip" k of X{p). A numerical evaluation shows that all poles 
are such that p' k < 0. At asymptotic times only those poles for which p' k = will be relevant. 
There is only one such pole, namely, po = —2iu leading to 



X(t) 



2 



Xi(0) -XiMo + /3 2 A xf 1 (0)o + xl(0) °Z l(2a,)o 



(50) 



Inserting this expression into Eq. ([45]), evaluating the integral, and returning to Eq. 
(|52]), we obtain the final result for the complex third-order nonlinear susceptibility (with r 
restored) : 



X fV) [xi(w)-Xi(3w)] 



Xi(0)o-XiMo + /3 2 A xi 11 (0)o+ Xl(0) ° 2 L X ; (2 " ) ° 



(51) 



Here Xi\ UJT ) is given by Eq. (p2|) with y = lot, and Xi^) by Eq. (p^). 

This expression may now be averaged over the distribution of relaxation times r as 
argued in the paragraph preceding Eq. (|39|) . In Fig. 4 the real and imaginary parts of XzW) 
are plotted as functions of temperature for several values of frequency. In analogy with the 
case of Xi( a; ) above, a VF behavior (|37D of r and a linear distribution of VF temperatures 
T has been used. The values of the parameters are again Jq/J = 0.9 and A/ J 2 = 0.001. 
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The real part Xi(T, ui) has a sharp peak near T ~ J, whose origin can be traced back to the 
function x^(co>) appearing in Eq. fl5T|). As in the linear susceptibility case, a strong frequency 



dispersion is evident. In the limit of small frequencies, i.e., ujtq <C 1, the behavior of x' 3 (T) 
is the same as in the static case studied by replica theory f?|. At high frequencies and low 
temperatures, x 3 (T,uj) may become negative due to the last factor in the numerator of Eq. 
(pT|), whose imaginary part changes sign. This effect cannot be observed easily, since this 
would require a measurement of the nonlinear susceptibility in the range where the absolute 
value of x' 3 (T, uj) is extremely small compared to its peak value. 

It is easily verified that in the limit u — > 0, Eq. (^1|) reduces to the static result (|27[) . 

It has been shown in the static theory |7J that a crucial quantity, which can discrimi- 
nate between the dipolar glass-like and ferroelectric behavior, is not Xz{T) but rather the 
rescaled static nonlinear susceptibility a 3 = Xz/Xi- m s pi n glasses without random fields 
a 3 (T) diverges at Tf, and in a relaxor it develops a peak near Tf, whereas in a ferroelectric 
with long range order a 3 decreases with decreasing temperature on approaching the critical 



temperature T c [O]. It has been suggested that a 3 (T) could also be extracted from the dy- 



namic linear and nonlinear susceptibilities by considering the following generalized function 
IJ,|24|: 



«< r -»>=;dSr (52) 

In Fig. 5, a' 3 (T, uj) is plotted as a function of temperature for the same set of parameters as 
in Fig. 4. Each of the three factors in Eq. (|52"D has been averaged over the linear distribution 
of To. Obviously, a 3 (T,u) develops a peak near Tf ~ J at all frequencies shown. On 
the high-T side of the peak, a' 3 (T) is independent of frequency and agrees with the static 
a 3 . Near the peak and on its low-T side, however, strong frequency dispersion appears. 
Recent experiments in PMN and PLZT || indicate that the high-T quasistatic part of 
a 3 (T) exhibits a crossover between the paraelectric-like decreasing behavior and a glass-like 
increasing behavior on approaching Tf. This type of behavior is characteristic of relaxor 
ferroelectrics. The crossover can be qualitatively described by the present dynamic SRBRF 
model, as shown in the inset of Fig. 5. The model also correctly predicts the onset of strong 
frequency dispersion in a' 3 (T, u) at low temperatures, which has been observed in both PMN 
and PLZT [§. 



VI. CONCLUSIONS 

We have presented a dynamic model of uniaxial relaxor ferroelectrics based on the re- 
cently developed static spherical random-bond — random-field (SRBRF) model 0. Following 



the theory of spherical models of spin glasses fl2| , |13|l the order parameter field is assumed 
to obey the Langevin equation of motion written in the representation of eigenstates of the 
random interaction matrix with the spherical condition being enforced at all times. The 
equations of motion are exactly solvable in the asymptotic limit where the relaxor system 
reaches an equilibrium state. The linear and the third-order nonlinear dynamic response 
functions have been derived. In the static uj —>■ limit these results are precisely equivalent 
to the static linear and nonlinear susceptibilities Xi an d Xz-, respectively, obtained earlier by 
the replica method [[7|. 
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In analogy with the static case, the dynamic theory does not predict a sharp transition 
into a dipolar glass-like state. Rather, in the case of weak random fields the third-order 
susceptibility shows a narrow peak at a temperature Tf, which mimics the freezing transition. 
Within the context of a dynamic model the freezing transition would correspond to the 
divergence of the longest relaxation time in the system. However, the dynamic SRBRF 
model contains no information on the behavior of the relaxation time r appearing in the 
equations of motion, and does not lead to the divergence of the effective relaxation time on 
approaching Tf. In order to describe the observed freezing transition one should therefore 
introduce a divergent behavior of r. This can be done, for example, by assuming a Vogel- 
Fulcher (VF) law for the temperature behavior of r in accordance with empirical findings 
HJfJ, however, this will suppress the dynamic response at all temperatures lower than the 
VF temperature To. We have shown that by introducing a probability distribution of VF 
temperatures T one can obtain linear and nonlinear response functions which remain finite 
at all temperatures, in qualitative agreement with experiments. The largest value of T has 
been set equal to the static "freezing" temperature Tf, which is determined by the random 
bond strength parameter J. 

The actual shape of Xi(T, u) and x 3 (T, to), where the bar means an average over r or 
equivalently over T , strongly depends on the probability distribution of relaxation times 
#(mr). 

Within the framework of the dynamic SRBRF model we have also calculated the scaled 
dynamic nonlinear susceptibility a' 3 {T, to) = x'^) /Xi(3u))xi(u) 3 , which allows one to dis- 
criminate between the ferroelectric-like and glass-like behavior of relaxors []5],|24|. In the 
quasistatic regime above Tf, a 3 (T,io) is practically independent of to and its temperature 
dependence shows a crossover between paraelectric-like and glass-like behavior on approach- 
ing Tf from above. This crossover behavior has recently been observed both in PMN and 
PLZT ||. The calculated shape of Xi{T, to) and XsiT, to), and hence of a 3 {T, to), strongly de- 
pends on the probability distribution of relaxation times g(lnr). In the present work we did 
not attempt to investigate in detail the effects of g(lnr) on the behavior of these quantities; 
however, we have shown that if one assumes a linear distribution of VF temperatures T , the 
predicted behavior of Xi(T,to) and x 3 (T, u;) is in qualitative agreement with experiments in 
PMN [gjniH and PLZT On the other hand, a 3 (T,tu) obtained in this manner has 

a peak near Tf ~ J and shows a strong frequency dispersion below Tf. The predicted peak 
in a' 3 (T, to) has not been observed experimentally |24]], suggesting that one should perhaps 
search for a more realistic distribution g(lnr). This problem will be dealt with in a future 
publication. 



ACKNOWLEDGMENTS 

This work was supported by the Ministry for Science and Technology of the Republic of 
Slovenia. One of the authors (R. P.) is grateful to Silvio Franz and to Boris E. Vugmeister 
for helpful suggestions. 



13 



REFERENCES 



[1] L. E. Cross, Ferroelectrics 76, 241 (1987); 151, 305 (1994). 

[2] B. E. Vugmeister and H. Rabitz, Phys. Rev. B 57, 7581 (1998). 

[3] R. Blinc, J. Dolinsek, A. Gregorovic, B. Zalar, C. Filipic, Z. Kutnjak, A. Levstik, and 

R. Pirc, Phys. Rev. Lett. 83, 424 (1999). 
[4] R. Blinc, A. Gregorovic, B. Zalar, R. Pirc, and S.G. Lushnikov, Phys. Rev. B 61, 253 

(2000). 

[5] V. Bobnar, Z. Kutnjak, R. Pirc, R. Blinc, and A. Levstik, Phys. Rev. Lett. 84, 5892 
(2000). 

[6] B. E. Vugmeister and H. Rabitz, Phys. Rev. B 61, 14448 (2000). 
[7] R. Pirc and R. Blinc, Phys. Rev. B 60, 13470 (1999). 

[8] V. Bobnar, Z. Kutnjak, R. Pirc, and A. Levstik, Phys. Rev. B 60, 6420, (1999). 

[9] V. Westphal, W. Kleemann, and M. D. Glinchuk, Phys. Rev. Lett. 68, 847 (1992). 
[10] R. Pirc, B.Tadic, and R. Blinc, Phys. Rev. B 36, 8607 (1987). 
[11] N. Korner, Ch. Pfammatter, and R. Kind, Phys. Rev. Lett. 70, 1283 (1993). 
[12] L.F. Cugliandolo and D. S. Dean, J. Phys. A: Math. Gen. 28, 4213 (1995). 
[13] L.F. Cugliandolo and D. S. Dean, J. Phys. A: Math. Gen. 28, L453 (1995). 
[14] J. Kosterlitz, D. Thouless, and R.C. Jones, Phys. Rev. Lett. 36, 1217 (1976). 
[15] D.C. Mattis, Spin Glasses (Springer- Verlag, Berlin, 1985), Chap. 2.12. 
[16] Y. Ishibashi, J. Korean Phys. Soc. 32, S407 (1998). 

[17] G.A. Samara, Ferroelectrics 117, 347 (1991); Phys. Rev. Lett. 77, 314 (1996). 

[18] Z. Kutnjak, R. Pirc, A. Levstik, I. Levstik, C. Filipic, R. Blinc, and R. Kind, Phys. 

Rev. B 50, 12421 (1994). 
[19] A. Levstik, Z. Kutnjak, C. Filipic, and R. Pirc, Phys. Rev. B 57, 11204 (1998). 
[20] Z. Kutnjak, C. Filipic, R. Pirc, A. Levstik, R. Farhi, and M. El Marssi, Phys. Rev. B 

59, 294 (1999). 

[21] H. Sompolinsky and A. Zippelius, Phys. Rev. B 25, 6860 (1982). 
[22] B.Tadic, R. Pirc, and R. Blinc, Z. Phys. B 75, 249 (1989). 
[23] R. Pirc (unpublished). 

[24] A.E. Glazounov, Ph. D. Thesis No. 1665, EPFL, Lausanne (1997); A.E. Glazounov and 
A.K. Tagantsev, Phys. Rev. Lett. 85, 2192 (2000). 



14 



FIGURES 



FIG. 1. Solution of Eq. (23) for the Lagrange multiplier z(T) obtained numerically in zero field 
(E = 0) and for A/J 2 = 0.001. Inset: Proof that z - (3 J > at all temperatures. 

FIG. 2. Real and imaginary parts of the linear susceptibility in the case of a single Vogel-Fulcher 
(VF) type relaxation time [Eq. (37)] as functions of temperature for several values of frequency, as 
indicated. Note that the response is strictly zero below the VF temperature To = J. 

FIG. 3. Real and imaginary parts of the linear susceptibility averaged over a linear probability 
distribution of VF temperature To, with < To < J. 

FIG. 4. Calculated temperature dependence of the real and imaginary parts of the third-order 
nonlinear susceptibility averaged over a linear distribution of To , using the same parameter values 
as in Figs. 2 and 3. 

FIG. 5. Temperature dependence of the scaled nonlinear susceptibility 
a' 3 (T,Lo) = X3(^)/x'i(3w)x'i( u; ) 3 obtained with the same parameter values as in Figs. 2-4. The 
inset shows the crossover from the decreasing paraelectric-like to increasing glass-like behavior in 
the quasistatic regime above the peak. 
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